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SUMMARY 


The  purpose  of  this  note  is  to  study  the  statistical  problems 
of  estimation  and  confidence  for  both  the  location  and  speed  and 
heading  of  a  target  following  a  linear  course  when  tracked  by  radar 
from  an  Observatory  Satellite  (0.  S . ) . 

Section  I  of  this  report  specifies  in  detail  the  assumptions 
upon  which  the  analysis  of  the  subsequent  sections  is  based.  The 
main  mathematical  results  are  contained  in  Sections  II  and  III. 

The  join’  distribution  of  the  estimates  for  target  position  is  derived 
in  Section  II  and  in  Section  III  the  maximum  likelihood  estimate  of 
targ  t  position  is  obtained  and  a  proof  given  that  it  is  normally 
d^  tributed. 

Section  IV  contains  a  detailed  comparison  of  the  mathematical 
models  of  this  and  the  earlier  report,  and  the  important  theorems  and 
formulas  relating  to  each.  Sections  V  and  VI  present  special  cases 
in  which  the  results  can  be  expressed  In  a  form  which  is  simpler  than 
the  general  results  of  Sections  II  and  III. 

The  last  two  sections  deal  respectively  with  confidence  regions 
for  target  speed  and  heading  and  the  possibilities  for  tracking  a 
target  in  a  dense  environment  given  certain  information  on  the 
distribution  of  target  speeds  and  directions. 

The  notation  used  in  this  report  is  such  that  much  of  the  detailed 
analysis  of  the  previous  work  Reference  [2]  applies  directly  with  at 


j 


most  a  reinterpretation  of  the  symbols  involved.  Those  results 
from  the  earlier  report  which  do  not  apply  directly  to  this  model 
are  proved  to  hold  for  this  case. 


I.  INTRODUCTION 


This  paper  is  a  continuation  of  an  earlier  report,  Reference  (2j, 
in  which  the  same  problem  for  the  location  of  a  target  was  considered 
under  slightly  different  assumptions  which  we  call  Model  I.  Both 
reports  make  similar  assumptions  regarding  the  errors  arising  from 
the  surveillance  of  a  target  by  an  0.  S.  The  difference  lies  in  the 
statistical  errors  made  in  the  determination  of  the  0.  S.  position 
from  the  ground. 

In  the  first  model,  the  position  of  the  0.  S.  was  determined 
independently  for  each  observation  of  the  target.  This  would  be  the 
case  when  the  observations  were  made  from  different  0.  S.'s  or  when 
observations  from  a  single  0.  S.  were  so  spaced  in  time  that 
independent  assessments  of  its  position  were  available  for  each  target 
observation. 

In  the  model  presented  here,  a  determination  of  the  orbit  of  the 
0.  S.  is  made  and  its  position  at  a  given  time  estimated.  Positions 
at  a  later  time  in  the  same  orbit  are  then  calculated  from  this  estimate. 
This  would  apply  to  situations  where  several  observations  are  made  on 
a  single  overflight  and  no  new  estimate  of  the  0.  S.  position  is  made 
between  the  first  and  last  observation.  Thus,  successive  estimates  of 
target  position  via  0.  S.  observations  are  independent  in  the  previous 
Model  I  while  they  art  dependent  in  Model  II  presented  here. 
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As  in  the  first  report  on  this  subject,  we  need  to  specify 
precisely  the  assumptions  on  which  the  analysis  is  based.  These 
are : 

1°  The  0.  S.  moves  at  a  known  constant  altitude  linearly 
over  the  plane  with  a  known  constant  velocity. 

2°  The  target  moves  in  a  straight  line  at  a  constant  but 
unknown  velocity  in  the  plane. 

3°  The  estimated  position  coordinates  of  the  0.  S.  over 
the  plane  as  determined  from  the  ground  at  a  given 
time  are  bivariate  normal  random  variables  with  known 
covariance . 

4°  The  estimated  position  coordinates  of  the  target  on 
the  plane  as  determined  relative  to  the  true  position 
of  the  0.  S.  are  bivariate  normal  random  variables  and 
successive  observations  of  such  relative  positions  are 
independent . 

5°  Time  between  successive  observations  can  be  measured 

with  sufficient  accuracy  so  that  errors  of  position  due 
to  time  inaccuracy  are  negligible. 

6°  The  parameters  of  the  covariance  matrix  of  the  observations 
of  target  position  relative  to  the  true  0.  S.  position  can 
be  determined  from  bearing  angle  and  range  data. 

Thus,  it  is  clear  that  Model  II  is  intended  to  supplement  Model  I 


and  not  to  replace  it. 


3 


In  general,  much  of  1  he  detailed  analysis  of  specific  cases 
applies  directly  to  both  models  since  we  have  made  a  special  effort 
to  utilize  the  same  basic  notations  in  both  reports.  Also,  in 
general  we  try  to  make  the  same  notational  conventions  of  Reference  [2] 
apply  here.  Some  of  these  are: 

(i)  Summations  will  only  be  over  the  indices  i,j,k,£ 
and  extend  from  1  to  n. 

(ii)  Creek  letters,  upper  and  lower  case,  will  in  general 
denote  respectively  random  or  non-random  points  in 
the  plane  as  vectors. 

(iii)  Upper  case  Roman  letters  with  or  without  affixes  such 
as  A,B, . . .  will  be  matrices. 

(iv)  Random  variables  will  be  indicated  by  upper  case  letters 
and  by  lower  case  letters  with  circumflexes. 

(v)  Numbered  equations  of  the  previous  Model  I  in  Reference 
[2]  when  referred  to  will  be  prefixed  by  a  I. 

II.  OBSERVATIONS  OF  TARGET  POSITION 

On  a  single  overflight,  an  0.  S.  may  make  several  observations  of 
the  position  of  a  target  and  our  first  problem  is  to  determine  the 
joint  distribution  of  the  observations  of  the  target  position  using  the 
fact  that  the  target  position  relative  to  the  0.  S.  is  subject  to 
observational  error  as  is  the  estimate  of  the  0.  S.  positions  relative 
to  the  ground.  Thus,  at  a  given  time  t^  we  assume  the  0.  S.  is  at 
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some  position,  say  p^,  and  the  target  at  some  position  with 

respect  to  a  given  coordinate  system.  However,  the  position  of  the 
target  as  observed  by  radar  from  the  0.  S.  is  subject  to  chance  error 
and  hence,  the  radar  estimate  of  the  target  position  from  the  0.  S. 
position  is  a  random  variable,  say  E^.  Now  a  radar  measurement 
from  the  ground  at  time  tQ  of  the  0.  S.  position  on  the  given 

coordinate  system  is  also  a  random  variable,  call  it  A. 

From  assumptions  1°  and  2°,  we  have  that  the  0.  S.  follows  a 
linear  path  in  the  plane,  say, 

P  *  i  +  et 

as  does  the  carget 

U  *  a  +  Bt . 

Without  loss  of  generality,  we  can  select  our  coordinates  so 
that  the  first  coordinate  of  e  is  in  the  direction  of  0.  S.  travel 
and  hence,  the  second  coordinate  of  c  is  zero.  Moreover,  by 
assumption  1°,  the  first  coordinate  of  i  is  known. 

Again,  following  the  general  mathematical  assumptions  3°  and  4° 
and  the  notation  of  our  previous  model,  we  have  that  H^,  for 
i  *  l,...,n,  and  A  are  bivariate  normal.  More  specifically,  E^ 
is  a  normal  variate  with  mean  and  known  covariance  matrix 

C^.  That  is,  E^  is  N(p^  -  p  ^ ,  C  and  thus  has  density  given  by 

exp  |  -*j(  S  -  Pi  +  _  +  V+|  £  e  R2. 


2n 
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Also,  A  has  mean  pQ  and  covariance  matrix  D  and  thus 


has  density 


jjp  exp  J  -*s(A  -  pq)D(A  -  Pq)  V  Ac  R2. 


The  data  that  we  obtain  from  the  0.  S.  yields  the  observations 
H  *  5^  +  A  for  i  ■  l,...,n  and  we  seek  the  joint  density  of 
H  *  ),  call  it  f.  Now  the  joint  density  of 

A)  is 


n  HH exp  -Hai  -  vi +  °i)ci<h  -  +  pi,t| 


f‘  JT  ‘'xpi  JS<»-P0>D(»-<>0)t 


from  which  the  density  f  may  be  found  by  using 


ni  *  £i  +  n  ■  <ni . V- 


Thus, 

f(n) 


|c  | 

/  ^  cxp  1  ~VUi  * x '  “i +  Di)ciui  ■ x  ■  ui +  pi,+ 


•  exp  |  -*s(»  -  eQ) D( A  -  en)  (  dX 


2  7T 


0 


where  the  integration  extends  over  the  A-plane.  If  we  denote 

u  *  x  •  po-  ■  ni  •  ui  +  pi  •  co  then 


f(n> 


/A 


U  l. 

ilnl^ 


i^riDi 


i-l  <  2tt  ) 


n+1  exP  j  +  J  >  dw  . 
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Now  in  the  exponent  we  have 


[•••]  ■  +  E(idC^u)  -  ~  ^  i  C  i  C  ^ } 

-  w(D  +  EC^o/  -  E^C^  wt  - 

Since  for  any  vector  t  and  any  matrix  A  we  have 

"t"  t"  "t*  *f* 

(ti)-T )  A(oj-t  )  -  iuku'  -  tAou  -  wAt  +  tAt  , 


if  we  let 


(2.1) 

then 

and 


-1 


A  -  D  +  EC  ,  t  -  D^C  A 


(...)  -  (u,  -  t)A(w  -  t)  -  tAt  + 


f  (n) 


1D!^  ^?lci ,&a 


(2tt)' 


exp  J  -M-tAt  +  Z^ici^1J 


/ 


lAj. 

2  TT 


exp[~4(w  -  t )A(oj  -  t)  ]du). 


Writing  c  -  (r  c  ),  we  now  want  to  find  a  matrix  E  so 

i  n 

that 

.i- 

(,E;  ■  ' ' r, ^ C i r. i  -  tAt  . 

Substituting  (2.1)  and  expanding  both  sides,  we  have 


2  V:ij'j  =  2 
i  j  i 


Wi  -JWVj* 
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using  the  fact  that  both  D  and  C,  hence  A,  are  symmetric 
Thus,  we  have 

/  _  i 

if  i  +  j 


(2.2) 


ij 


-v'S 


ci  -  V'lci 


if  i  -  j, 


from  which  one  checks  that  E  *  ^ j  i  anc*  hence  that  E  is 

symmetric,  itself.  Denoting  +  pq  then  r7 ^  ~ 

and  C  *  n  -  u,  we  have  the  fundamental 

T'.EOREM  1:  The  distribution  of  the  observations  of  target  position 

H  =  (1L,...,H  ),  obtained  from  the  0.  S.  is 
1  n 


(2.3) 
where 

(2.4) 


f  (n) 


(20 


exp  [ - *5 ( r,  -  u)E(n  -  u)  ]  n  c  R 


2n 


MIJIO 

E|  «  -  .  Thus,  we  conclude  that  if  is 


distributed  as  N(u,E-^).  That  is,  H  is  normally  distributed  with 
mean  u  and  covariance  matrix  E 


To  complete  the  proof  of  Theorem  1,  we  should  verify  Relation  (2.4) 


From  (2.2)  we  have 
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E  - 


-1 

1 

1 

„  .-1 

C .  A  C, 

-  C  A 

c. 

-  c_A 

C,  ... 

-  C  A 

1  1 

2 

1 

3 

1 

n 

-1 

1 

<M 

1 

,  -1 

C.A  C. 

-  c„a~ 

C' 

-  C„A 

c„  ... 

-  C  A 

1  2 

2  2 

2 

3 

2 

n 

-1 

1 

1 

-1 

C  A  C., 

-  c2a" 

c. 

C..  -  C  A 

C„  ... 

-  C  A 

1  3 

• 

• 

• 

• 

• 

• 

3 

3  3 

• 

• 

• 

3 

• 

• 

• 

n 

• 

• 

• 

C,A_1C 
1  n 


-  C0A  C 
2  n 


C„A-1C  ...  C 


C  A  XC 
n  n 


^  |  ^ 

Now  we  may  factor  from  the  i  column  for  each  i  *  l,...,n  and 

subtract  the  first  column  from  the  2nc^,  3rc*,...,  last  column  to  obtain 


Ei-n  ict 

i  =  1 


I  -  A_1C, 


-  A_1C, 


-  A_iC. 


-I 

I 

0 


-I  ...  -I 

0  .  .  .  0 

I  .  .  .  o 


-  A_1c 


0 


0  .  .  .  I 


th 


Multiplying  the  i  column  by  A  and  adding  to  the  first  column, 


we  get 


ni 


1  -  ^A_1C 


0 

0 


-I 

I 

0 


0 


0 


-I  ...  -I 

0  .  .  .  0 

I  .  .  .  0 

•  •  «  •  • 

•  ■  •  •  • 

•  •  •  •  • 

0  .  .  .  I 


I  -  EA-1C 
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To  complete  our  check,  we  show  that 

|l  -  ZA~1Cil  «  |A_1A  -  IA_1C1|  -  | A*1 1  | A  -  ECj  -  jy-  . 

III.  THE  MAXIMUM  LIKELIHOOD  ESTIMATES  AND  THEIR  DISTRIBUTION 

Theorem  1  tells  us  that  a  single  observation  of  0.  S.  position 
from  a  bivariate  normal  distribution  relative  to  the  ground  combined 
with  n  observations  of  target  position  which  are  normally 
distributed  relative  to  the  0.  S.  position  yields  a  joint  normal 
distribution  of  observed  target  positions  relative  to  the  ground. 

Our  next  problem  is  to  find  a  best  estimate  of  target  course  using 
this  normal  distribution  of  error  and  the  assumption  that  the  target 
is  moving  linearly  at  constant  velocity. 

Thus,  we  want  to  find  the  maximum  likelihood  estimate  (m.l.e.) 
of  future  target  positions  call  the  estimate  Now  by  the 

known  properties  of  the  maximum  likelihood  estimates  we  know  that 

Ut  ”  a  +  Bt, 

where  a  and  8  are  the  m.l.e.  of  a  and  8  and  to  complete  the 
picture  we  need  to  know  the  distribution  of  (a,B). 

Now  we  see  that 

Ui  =  Ui  ~  Ci  +  P0  *  +  ct0>  +  “  c)ti. 
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If  we  define 

Y  *  a  +  etQ,  x  ■  6  -  e 

we  have  a  model  similar  to  that  analysed  previously  and  if  we 
obtain  the  maximum  likelihood  estimates  of  y  and  k  ,  then 
we  can  easily  obtain  the  m.l.e.  of  a, 6. 

The  likelihood  function  is 

L  -  Ln  |  E  |  -  nLn(2n)  -  2  ~  ui*Eil*nl  "  '• 

ij  J  J  J 

Let  F  ■  with  and  let  us  for  this  argument 

introduce  the  notation  for  the  coordinates  of 

Y  ■  (Yj.Yj)*  *  * 

then 
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holds  for  all  j  *  l,...,n  with  6  the  Kronecker  delta,  we 
obtain 


(3.1) 


3y  "  22Wk  +  6kEijCj) 
k  iJ 


k  -  1,2 


(3.2) 


jL  _  1 
3k 


k  ij 


2  2  <tJ^iElj6k  +  ti<5kEij;j)‘ 


Thus,  we  have  from  (3.1),  by  realizing  that  the  second  term 
is  a  scalar  and  thus  equal  to  its  own  transpose, 

2(ii +  VEij4k  ■ 0  k’1-2- 

u 

or  writing  these  two  equations  in  matrix  form 

+  Cj)Eij  <6l*62)  *  (0‘0)’ 

ij 

But  since  (*1*^2)  =  (3.1)  reduces  to 

(3.3)  2^!  ♦  t,)E  =  (0,0). 

ij 

Now  from  (3.2),  following  similar  steps,  we  have 

2<Vi +  Vj’VJ  ’ 0  k ' 1,2 

ij 

from  which  we  obtain 

2(tj;i  +  tiCj)Eij 

xj 


(3.4) 


(0,0). 
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But  wo  no to 


i  +  'j  ’  "i  '  '  '  -ti  +  "j  ‘ '  ‘  ,lj ’ "t  +  r,j  ‘  2<  ‘ '  (ti  +  V 


Y<  +  Y'J  '  ‘j"!  •  V‘j  •  ‘Vj  +  Y'J  ‘  V‘i  '  ‘Vj 


V‘t  ♦  t|nj  -  ,<t,  ♦  tj>  -  2«Vj. 


Hence  by  substitution,  wo  finally  obtain  from  (3.1)  the  equation 


(3.5) 


y  (r,  +  r,  )K 

i  j  ij 


2t  2  Kij  +  *  ^  (b  +  V: 


ij 


and  f  rom  (3.2) 


(3.6) 


2  <t.v.  +  t  .r  ,)K.  .  =  >  2  <t.  +  t  .)E.  +  2*  V  t  t  F. , 

Jj  J  1  i  J  i J  ff  i  J  ij  fj  i  J  ij 


Writing  those  in  matrix  notation,  wo  have 


(3.7) 


(t,>  ) 


‘  1 1 


r 

21 


‘12 

22 


(’ 1*^2'  * 


V  =  V  y 

—  -fa-  i  i 

11  ii  J 


2  -2 

12  21 


V 

— 

i  1 


t  .  +  t  . 

_L _ L 

> 


j 


where 


Now  it  we  define  the  matrix  F 


(V 


by 


F  =  hilh. i 

ij  2 


E  +  F 

wbicli  we  write  in  matrix  form  F  =  — - -  with  the 

definition 


F.°  =  (E°  )  where  E°  =  E.  .  , 
ij  ij  Ji 


then  we  can  see,  by  rearranging  some  summations. 


(3.8.1)  2  *  2  Fn 

11  ij  J 


(3.8.2) 


V 

4m* 

ij 


t 


i 


F.  . 

ij 


(3.8.3) 


2  “  2  t  .  t  .  F .  . 

.  *■!  l  j  ij 

22  l  j 


and  in  the  same  manner 


(3.8.4) 


2 

i  j 


-  .  F  .  . 
i  n 


(  1.8.3) 


t  , 


.  +  t  . 
i _ ! _ 1 

2 


obvious 


Now  if  we  denote 
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( 3.9.1) 

(3.9.2) 


wr l t ten 

( Y  « * 

•lSll  + 

’  2S2  1 

c  s  *f 

1  S 

112 

*2  22 

S21  *  S 

12* 

1  ’ v  2 


( C )  *  ( n n  ) 
I  /  L  n 


Y  t  F 
j  J  1J 


V  t  F 

7  J  nJ 


\ 

/ 


=  nU 


with  ri,U  do  fined  in  t  he  obvious  manner.  Since  by  Theorem  1, 

'  is  N(v,K  S  and  linear  transformations  of  normal  variables 
are  again  normal,  it  is  clear  that  .  =  is  bivariate 

normal.  To  be  precise,  by  P.  19  of  Reference  [1],  we  have  that  i* 
is  N(vU,U  E  V).  Since  ( C , *  )  =  .S,  we  know  by  a  similar 
argument  that  (0,>)  is  N(vUS,Q)  where 


(3.9.3) 


()  =  s^u 1 1;“ 1  is . 


We  will  now  prove  that 


(  3  .  10)  vl'S  =  (  v,*  )  • 

To  do  this,  we  make  use  of  the  identities  which  follow  from 


o.  lo.oi) 


is 


‘  1  i S 1 1  +  12S21  “  If  '  1 1 S 1 2  +  12S22  =  C' 


A21S11  +  ‘ ' 2 2 S 2 1  *  2 1 S 1 2  +  '22S22  *  1‘ 


The  first  component  of  the  vector  vlS  is 


2<v  +  '■‘l>r.jsnt2  <v  +  •S>tJFiJS2i 


ij 


ij 


Y°  HS11  +  l\2S2>'  +  ’  (I21S11  +  ‘  2 2 S 2 1 }  “  '* 


Tlie  second  component  of  the  vector  vt’S  is 


~  <v  +  •ti)Fijsi:  +  2  (v  +  •ti;tjFiJS22 


Y (' 1 1S12  +  ‘  12^22  ^  +  *  ^  ?l Sl ?  +  ‘ 


2112  22  22 


Tli i s  proves  (3.10).  Now  under  the  additional  assumption  that 


(3. 10.1) 


E°K  [le  =  E 


we  can  show  that 


(3.10. 2) 


0  =  S. 


If  we  write 


F  f 

11  In 


U  =  FT0  = 


n  1 


\  /'•  T 


nn 


/ 


\  '• 


t  I 
n 
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where 


(3. 10.3) 
then 


TJ  '  <,'V) 


U+E_1U  -  TFE_1FT° 


T  -  (Tr...fTn) 


since  F  *  F  we  have 


4FE**1  F  -  (E  +  E°)E-1(r.  +  E°) 


=  E  +  E°  +  E°  +  E°E-1E° 


2F. 


Now 


TFT0 


CjFij 


t  .  t  F 

i  J  ij 


S 


-1 


and  the  assertion  follows  since  S  =  S. 

We  have  just  obtained 

THEOREM  2:  The  maximum  likelihood  estimates  (a,B)  defined 
by  (x  =  y  -  t  t  ,  P  =  <  +  t 

where  (9,<)  are  defined  by  (3.9.1),  (3.9.2)  and  are  N((a,6),Q) 
with  Q  given  in  (3.9.3). 

From  this  we  have  the  immediate 

COROLLARY  3:  The  maximum  likelihood  estimate  of  the  true  target 
position  u  at  time  t  is 
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=  (i  +  r  t 

wliich  Is  distributed  as  ) 

where 

(3.11)  Bt  «  Qu  +  t  ( Q  j  2  +  Q21)  +  t2Q22- 


We  will  now  prove 


THEOREM  3:  The  estimate  f;  and  its  covariance  matrix  are 

invariant  under  location  and  scale  change  in  time. 


PROOF:  Let  t'  =  at  +  b  be  a  linear  transformation  of  time  and 
use  primes  affixed  to  matrices  to  denote  that  matrix  computed  at 
the  new  times 


(3.12) 


I  _ 


11 


ir 


i 

"21 


=  a 


21 


+  b).' 


ir 


•  » 

"22 


2.. 

=  a  i. 


22 


+  2ab2 


2. 


21 


11* 


Now 

Ot  =  n  +  Bt*Y  +  iCt*r(t”  tQ) 

where  i  represents  the  velocity  of  the  0.  S.  per  un  i  _  i  :  l.m  . 
is  thus  sufficient  to  show  that  y  +  >“t  is  invariant.  In  our  previous 
notat ion 


(YV)  =  r  F  T 0  S 


where  r  F  is  independent  of  time.  Hence 


S11  +  tlS21  +  tS12  +  ltlS22 


Y  +  ►“  t  =  F 


b , ,  +  t  S  +  tS  +  tt  S 
11  n  21  12  n  22 


It  is  sufficient  to  show  that  foi  each  j  *  1,. 


•  • 


,n  that 


(3.U)  IjS1i  "  Sn  +  ljS21  +  liSl2  +  titjS22 


is  invariant. 


From  the  equations  (3.10.01),  (3.10.02),  we  have 


„-l  ,  ..  -1 

22  "  22  “  ^21  11  12 


S12  "  ”^Hri2S22 


Q  m  Q 
21  12 


Su  ■  ’ll*1  -  h2s21> 


There  follows  by  substitution  from  (3.1.2) 


c,-l  2,,.  _  -U  .  2-1 

S22  3  (122  "  '21  11^12  3  S22 


hence 


S^2  =  -y  S22*  Likewise,  we  find 
a 


*.  _  Ill  b  ,  _  I  ,  _L  c 

12  a  a  22*  21  ~  a  21  2  22 

a 


Sll'Sll-I(S21+S12)+i!2S22 

a 


Now  we  calculate 


T?'Tf  ■  S11  -  !<S21  +  S12>  +  s22  +  (atj  +  b)(7  S21  -  - 

J  a  J  a 


+  (at.  +  b)  (~  S12  -  -y  S22^  +  ^ati  +  b^at^  +  b) 

a 


expanding  and  simplifying  siiows 


S11  f  ljS21  +  CiS12  +  titjS22  =  1 J S 1 i 


Thus,  we  have  shown  that  a  +  «'  t  is  invariant  under  sc  ale  and 
location  change  in  time. 


To  see  that  B  is  also  invariant  under  the  same  transformation 
t 

of  time  we  merely  note  that  if  we  write  T  *  (l,t^I)  then  by  making 
use  of  the  notation  of  (3.10.3)  and  the  fact  that  T  =  FT0  we  have 


Bt()  =  ToOTo 


(T()ST  )  F  K 


-1 


F(TST()) 


but  T^ST  is  itself  invariant  by  our  previous  argument. 


IV.  A  COMPARISON  OF  MODELS  I  AND  II 

We  have  proceeded  thus  far  without  a  direct  comparison  with 
Model  I.  As  in  Section  II  of  Reference  [2),  we  assume  that  at  a 
given  time  t.  we  have  the  0.  S.  at  point  c 


0.  S.  direct  ion 
of  travel  at 
known  velocity 


FIGURE  1 
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observing,  at  bearing  angle  ami  range  R^,  the  target  position 

j  .  The  observed  position  of  the  target  Is  the  random  variable  •  ^ 
which  by  assumption  is  N(ij^  -  t  ^y  following  the  same 

argument  used  in  Section  II  of  Reference  [2),  we  have 


Vi 


(4.1) 


-1 


Vi 


where  exactly  as  (1.2.3)  we  have 


(4.2) 


Vi 


2  2 

o  -  ■ 
u  .  V  . 

cov(x^y^)  =  — - - -  sin  2“ 


•> 


i* 


Following  (1.7. 1.2)  and  (1.7.  1.3)  we  have 


(4.3) 


2  2  2  2.2 
o  -a  cos  v .  4-  o  sin  *• . 

X  .  U  .  1  V  .  1 

1  1  1 


2  2  2  2  2 
(4.4)  =  o“  sin"''.  +  o  cos4"1).. 

v .  u,  l  v .  l 

i  i  l 


Now  at  time  t.,  n  is  the  standard  deviation  of  the  range  error 
1  Ui 

and  o  is  the  standard  deviation  of  the  azimuth  error  and  both 
v . 
l 

are  known  functions  of  the  given  operational  range  R^.  It  is 
clear  that  we  also  need  the  assumption  6°  that  the  parameters  of 


can  be  determined  by  the  bearing  angle  6^  and  the  range  . 
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Now  we  also  assume  as  in  Model  I  thai  the  matrix  I)  is 
diagonal  and  known,  namely 


The  difference  in  the  models  is  essentially  in  the  different 

densities  of  the  set  of  observations  n  ■  (n,,...,ri  ). 

1  n 

In  Model  I  we  have 


n  S 

n  im  k 

f(n)  =  i_J - -  exp  |  -'-s  2  "  ui)Ai(ni  ‘  } 

(2tt)  1 

where  =  C^  +  D  *  for  i  =  l,...,n. 

On  the  other  hand,  the  density  of  n  for  Model  II  is  given  in 
Theorem  i,  Equation  (1.10),  which  may  be  compared  with  that  above. 


V.  SPECIAL  CASE  OF  THE  DISTRIBUTION  OF  RADAR  ERRORS  CONSTANT  IN  TIME 

For  this  case,  we  assume  that  neither  the  orientation  of  the  0.  S.' 
direction  of  travel  nor  its  position  relative  to  the  target  plays  a  role 
in  determining  the  covariance  matrix  of  the  observation;  thus,  we  have 


C 


i 


n . 


C 


for 


i  =  1 


I  •  '  '  » 
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Now  from  (2.1)  wo  have 


A  ■  I)  +  nC 


and  f  rom  (2.2) 


E 


ij 


6  C  -  CA_1C 


Then  Q  =  S  and 


E 


ij 


1:,1  "2  Em  “2  ‘  nCA  1(:)  =  nC  -  nZCA_AC, 

ij  ij  i 


2.,  -1, 


“12 


l2\  =  2  t  E  *2  (tic  '  nt  CA"1(:)  *  ntC  -  n2tCA  1C, 


22 


CA  *0, 


-  2  t  .  ( t . C  -  ntCA  *C)  =  nt^C  -  (ntJ^CA  AC, 
i  1  J 


2„  -1, 


.  =  2  r'iEM  =  2  (n.C  -  nnCA  *C)  =  nn(C  -  nCA  C) , 
ij  J  i  1 


__  _  i  _  a*  “1 

ii<2  =  2  n.t  F...  =2r'i<t-C  “  nt^A  c)  =  nr.tC  -  n  tCA  C 


ij  -  j  ij  i 


By  the  invariance  theorem,  we  can  choose  the  time  origin 
so  that  t  =  0.  Hence,  we  have 


2  J 


2  -1 

/.  1 1  *  n(  -  n  (  A  (  , 


(5.D. 8) 
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'  12 


22 


nt  , 


Now 


*  nY-  j  j ,  ,  2  =  nt'iC  , 


.-1 

’ll  =  'll’ 


s  =  o  =  s 

12  21  * 


-1 


>22  =  i  22’ 


(5.0.9) 


“t  '  Su  +  '\2  *  :U  +  lh'l2' 


and  we  see  from  the  above 


Y  =  ’1S11  =  '■'» 


=  S 

*2*12  —  * 

t 


(5.1) 


+  .  . 


These  two  equations  may  be  compared  with  equations  (6.4.1)  and 
(6.4.2)  of  Reference  [2)  with  t  =  0. 


From  (5.0.8)  we  have 


2U 


) .  *  nCA  *(A  -  nC)  -  nCA  ^1), 


;  -MbV.IcSd-1, 

11  n  n 


>  -1  -  -L.  c-i 
22  - 
nt 


Now  substituting  into  (5.0.9),  we  have 


(5.2) 


B  -  -  (1  +  —  )C  1  +  D  1 
t  n  — 

t 


Thus 


,  for  this  case  we  have  easy  formulas  for  the  estimate  0 


and 


the  covariance  matrix  B^. 

Of  course,  if  we  assume  further  that  C  is  diagonal,  since  D 
is,  the  calculation  of  B^  is  immediate. 

We  remark  that  a  comparison  may  also  be  made  between  the 
two  models  by  examining  (5.2)  above  and  (6.3)  of  Reference  [2], 
remembering  that  in  Reference  [2]  A^*C^+D^. 


VI.  THE  SPECIAL  CASK  OF  SYMMETRIC  OBSERVATIONS 

In  this  case  we  have  the  n  observations  symmetrically  spaced 
in  time.  Using  the  invariance  theorem,  we  can  choose  the  origin  so 
that  we  have  the  relations 


t.  +  t  .  *  0 
l  n+l-i 


•.  +  o  ,  . 
l  n+l-i 


l  ~  1 , . .  .  ,n 
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for  the  observation  times  and  the  bearing  angles.  Now  if  we  set 


a 


i 


then 


XiVi 


t 


C 


-1 

i 


Letting 


.-1 


3iCi 


we  can  write 


e 


i 


and  in  particular  we  have 


e  . 
l 


sin  2^  . 
i 


From  the  above  relations  we  see  that 


di  =  dn-i+l' 


fl  fn-i+l’  1  "  1 . . 


ei  =  “en-i+r 


so  that 


Ci  +  ('n-i+l 


r  °) 


2b 


!  Inis  wt*  see  that  the  matrix 


Di  ■  i(:i 


is  diagonal  as  in  A  *  D  +  Dj .  A  matrix  of  the  form 


P  = 


(I,  :■) 


we  shall  call  contra-diagonal.  Hence,  the  matrix 


B2  ■ 


is  contra-diagonal.  Since 


Eij  ■  V.  -  v'v 


we  have 


i:n  -  js  Eij  ■I'fuVf'A1 


-2(t,  - 


=  I)1  -  l)  A_1D  . 


Since  the  inverse  product  and  difference  of  diagonal  matrices  are 
diagonal,  we  have  that  is  diagonal.  Now 


21  =  n  tit,j 


j§  *.<V.  -V'v 


-  tiSA'lui)  '  u2  -  D2A'‘V 
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And  is  diagonal  since  the  product  of  contra-diagonal 

matrices  is  diagonal.  Calculation  from  this  point  proceeds  in 
a  straight-forward  manner.  Since  we  have 

S21  *  (Z21  '  E11  21I22)  ’ 
by  equations  (3.10.01) 

(6.4)  Sx  =  -z21Z22S21’  S22  =  ■Z21Z11S12’ 

we  can  obtain  the  structure  of  S.  Since  contra-diagonal 

so  is  its  inverse,  and  in  particular  if 


Since  and  1^2  are  diagonal,  ^21^  ^  *s  t*ie  difference 

two  contra-diagonal  matrices  and  hence  is  contra-diagonal.  There¬ 
fore,  we  conclude  that  S^  is  contra-diagonal.  Thus,  we  see  by 
(6.4)  that  S^  and  S^2  are  both  diagonal.  This,  of  course,  is 
essentially  what  we  have  computed  in  Section  7  of  Reference  [1], 

Note  that  in  this  case  in  general  we  do  not  have  symmetric, 

as  we  had  for  Model  I. 

Thus,  we  see  that  all  elements  of  S  are  zero  except  those  on 
the  diagonal  or  contra-diagonal.  This  knowledge  may  be  of  help  in  the 
computation  of  Q  for  in  general  for  this  case  we  do  not  have  that 


Q  =  S. 
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VII.  ESTIMATION  AND  CONFIDENCE  FOR  SPEED  AND  HEADING 

The  parameters  with  which  one  is  usually  concerned  when 
tracking  a  moving  ship  target  are  the  speed  and  heading,  and  in 
this  section  we  shall  consider  t lie  accuracy  with  which  they  can 
he  estimated.  In  both  Models  I  and  II  we  have  obtained  the 
estimate,  fi(t)  ■  i  +  bt,  which  we  now  write  in  coordinate  form 


u(t)  =  (mj ( t ) ,  m^ ( t ) ) . 

This  is  an  estimate  of  the  target  position  ai  any  time  t 

wht n  it  follows  a  linear  path  in  the  plane.  Now  the  true  velocity 
s  of  the  ship  is 

s  =  ( (mj (t) ] 2  +  jm^ ( t ) ) 2  T5 
and  the  true  heading  h  is  the  angle 


arctan 


ni^(l)  -  m2(0) 
nij  ( 1)  -  m1  (0) 


plus  a  constant  depending  upon  the  signs  of  the  numerator  and 
denominator  of  the  argument  of  the  arctan  function.  A  little  more 
detailed  analysis  shows  that,  since  H  =  (b^jb^), 


/r 

=  V  bj 


+  b 


2’ 


and 
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b2 

0 

if 

b. 

s  o, 

b2  >  0 

h  -  arctan 

h 

1 

TT 

if 

bi 

<  0 

I 

i  2t! 

if 

bi 

>  0, 

b2  <  0 

where  arctan  is  the  principal  value  of  the  inverse  tangent  and 
takes  values  on  ( — tt / 2  , tt / 2 ) . 

These  relationships  are  seen  to  be  simply 

(7.1)  b^  *  s  cos  h,  b2  *  s  sin  h. 

Thus,  by  analogy  we  have  the  equations 

(7.2)  b^  *  V  cos  <J>,  *  V  sin  $ 

defining  the  random  variables  V  (for  velocity)  and  $>  for  the 
heading  angle,  which  are  estimates  of  the  true  speed  s  and  the 
true  heading  h.  (Here  we  violate  our  convention  to  use  Greek 
letters  for  vectors  in  order  to  retain  the  more  universal  one  that 
Greek  letters  are  angles.) 

We  were  able  to  show  that  the  estimators  (5,8)  are 
N(a,8),Q)  (Theorem  1  here  and  Theorem  2  of  Ref.  [2]),  where 


is  known.  Now  from  known  results  for  normal  variates  (p.  24  of 
Reference  (1))  it  follows  that  8  is  N(8,Q22). 
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'Pit*  Joint  density  of  $,V  is  found  by  making  a  transformation 
to  polar  coordinates  in  the  probability  density  of  N(B,Q22)*  This 
is 

i  v  ,  ,  —  It,  0  <  v  <  ® 

(7.3)  g  ( t ,  v  |  h ,  s )  =  - —  exp  {  } 

2v  ^22^  0  '  <  2tt 

where  were  we  have 

;  *  (v  cos  <$>  -  s  cos  h,  v  sin  <}>  -  s  sin  h) . 

This  density  can  be  used  to  study  the  distribution  of  velocity  and 
heading  estimates  that  could  arise  under  special  circumstances  such 
as  infrequent  headings  and  high  velocity. 

One  might  desire  separate  confidence  intervals  on  the  heading 
and  on  the  velocity.  However,  if  we  proceed  to  find  the  marginal 
densities  of  v  and  ^  from  (7.3),  we  see  that  each  density  has 
both  parameters  h  and  s. 

Thus  a  confidence  interval  for  the  velocity  could  be  constructed 
only  if  we  knew  the  true  heading  h.  Likewise,  a  confidence  interval 
could  be  found  for  the  true  heading  if  we  knew  the  true  velocity. 

The  presence  of  the  nuisance  parameters  prevents  us  from  obtaining 
confidence  intervals  separately  when  both  parameters  are  unknown. 

However,  we  can  obtain  a  joint  confidence  region  for  (h,s), 
which  is  somewhat  inefficient,  as  follows.  From  well-known  results 
on  the  Chi-square  distribution  of  the  quadratic  form  of  normal  variates, 
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we  have 

(7.4)  P[(8  -  e)Q"J(B  -  8)"  1  x^P))  =  1  -  P 
2 

where  X£^P^  was  Riven  in  Equation  (5.2)  of  Reference  [2]. 

Now  if  one  draws  the  ellipse  centered  at  6,  letting  C  =  (x,y) 
here  be  an  arbitrary  point  in  the  plane,  we  can  specify  the  ellipse 
as 

(7.5)  W  =  U  -  R)q‘*(c  -  8)+  1  x^pH- 

This  defines  a  lOO(l-p)  per  cent  confidence  region  for  6,  that  is, 

P 1 6  e  W]  *  1  -  p. 

Now  let  us  construct  the  smallest  area  in  polar  coordinates  which 
is  the  Cartesian  product  of  intervals  and  contains  U.  Let  us  call 
it  W*. 


FIGURE  2 
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2  /  2  2 
If  c;*(x,y)eR,  we  define  mod(0  *  r  x  +  y  ,  and 


arg(0  * 

arctan(y/x).  Then  we 

set 

(7.5.1) 

4>  -  sup  arg(0  , 

$  «  inf  arg(0 

SeW 

;eW 

(7.5.2) 

V„  *  sup  mod(0 , 

V.  *  inf  mod(c) 

te  W 

L  CeW 

and  define,  in  turn 

(7.3.3)  W*  =  1(9, r)  :  *  <  6  <  <,2*V1<  r  < 

Now  (r  cos  0,  r  sin  0)  e  W  implies  (0,r) e  W*. 


V* 


Thus 


P[(h,s)  e  W*]  >_  P[(s  cos  h,  s  sin  h)  e  W]  *  P [ 6  e  W]  ■  1  -  p. 


The  task  to  which  we  now  address  ourselves  is  the  determination 
of  random  variables  ,  i  =  1,2. 

The  boundary  of  the  elliptical  region  W  given  in  (7.5)  is 

(7.6)  (x  -  bj,  y  -  b2)Q22(x  -  bj,  y  -  t>2)  ’  -  X^P)* 

If  we  write 


by  symmetry  s^2  *  s2i’  anc*  we  use  German  letters  to  denote 


the  constants  defined  by 
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a 


sll 

2  * 

x2(p) 


b 


12 

2  * 
X  2  ( P  > 


C 


S 

22 

xj(p) 


then  we  may  write  (7.6)  as 


(7.7)  <l(x  -  b^)2  +  2b(x  -  b^My  -  t>2)  +  C  (y  -  t>2)2  -  1. 


We  want  to  find  the  maximum  and  minimum  of  the  functions 
(7.7.1)  f(x,y)  =  y/x,  g(x,y)  x2  +  y2 


subject  to  the  restriction  (7.7). 


Although  one  could  proceed  by  analytic  geometry,  or  by  the 
method  of  Lagrange  multipliers,  they  both  lead  to  the  solution  of 
a  quart  ic  equation.  A  simpler  method  seems  to  be  the  following. 
Make  a  transformation  of  (7.7)  to  the  variables  (u,v)  via 


(7.7.2) 


where 


x  *  bj  +  u  cos  0q  -  v  sin  6^ 
y  *  b2  +  u  sin  0^  +  v  cos  0^ 


2b 


tan  20o  ■  C  -  0 


Titus,  we  want  to  extremize  the  function  of  (u,v)  obtained  by 
substituting  (7.7.2)  into  (7.7.1),  subject  to  the  restriction 

2  2 

d^u  +  =  1, 


where  as  in  Equations  (5.4.1),  (5.4.2)  of  Reference  [2]  we  have 


2  2 

*  fl  cos  0q  +  b  sin  26^  +  c  sin  0Q 

2  2 

m  0  sin  Oq  -  b sin  +  C  cos  0Q. 


Now  if  we  let  cos 


/  a  ^  u,  sin  «  /  a  2  v  and  set 


cos  0, 


sin  0, 


i  -  1,2 


our  problem  becomes  that  of  extremizing  the  functions 


(7.7.3) 


and 


f($) 


b^  cos  $  +  fl 2  sin  $  ♦  b2 
0^  cos  4>  -  ^2  sin  $  +  b^ 


(7.7.4)  g(<*>)  »  Ca  1  cos 


-  b2  Sin  4  +  b^)‘ 


+  (bx  cos 


+  a2  sin  <t>  +  62) 


for  the  range  of  values  0  <  <f>  <_  2ir.  This  can  be  done  by  elementary 
techniques  of  calculus.  One  can  show  in  a  straightforward  manner  that 


f  *  ( )  “0  if  and  only  if 

(7.8)  a  +  6^  sin  41  +  82  cos  41  *  0 

where 
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with  capital  letters  those  random  variables  which  are  functions 
of  the  estimates. 

Now  let  be  the  unique  random  angle  for  which 

sin  $0  ■  B^/B,  cos  ■  Bj/B 

where 

B  -  +  B*  . 

Then  (7.8)  becomes  simply,  by  using  the  trigonometric  identity 
for  cosine  of  a  sum, 

cos(#0  -  $)  «  -a1/B. 

Letting  0  be  the  angle  on  (0,u)  such  that  0  ■  arc  cos(-a^/8) 
we  have,  since  cosine  is  an  even  function. 


°1  "  *0  -  °’ 


*o*0 


as  the  two  solutions  of  (7.8).  Thus,  we  have  as  the  limits  on  the 
true  bearing  angle  as  in  (7.5.3),  the  angles 

♦  j  »  arctan  f(fi^)  $2  *  arctan  f(n2). 

This  accomplishes  the  location  of  local  extremum  of  f.  Graphical 
examination  of  the  function  on  (0,2tt)  may  be  necessary  to  see  if 
these  are  the  true  extremum. 
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Proceeding  with  the  extremizing  of  g  we  see  in  a  manner 
exactly  similar  for  f  that 

g'($)  *  0  if  and  only  if 


(7.9)  sin  2<fr  +  Cl^  cos  2$  +  Cj  sin  <p  +  cos  <t>  *  0 


where 


I..  .  2  2  tb  .2 

2a^  n  ^  "°2  +  ^1  ~  ^2 


a  =  a  b  -  b  a 

2  12  *12 


C1  ■  bl«1  +  b2kl 


C2  ’  blb2  "  62°2 


and  again  upper  case  letters  denote  random  variables.  Unfortunately 

the  solutions  of  (7.9),  call  them  r  ,T  ,  must  be  found  in  general 

1  2 

by  some  numerical  technique.  However,  this  is  easily  accomplished  by 
a  computer  so  that  we  proceed  no  further. 


Having  found  the  solutions  1"^  i  *  1,2  the  limits  on  the 
velocity  which  were  sought  can  be  written 


Vx  -  min(  /I(iy,  /iTrJT) 

(7  9.1) 

V2  -  max(  /gTrjT,  /gTry). 


In  the  special  case  that  b  -  0,  (7.9)  simplifies  somewhat. 

This  circumstance  occurs  whenever  the  observations  are  drawn 
symmetrically  on  each  overflight,  so  that  it  is  of  practical  Interest. 
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In  this  case  0^  ■  0.  *  0,  (l^  ■  C 


a 


1 


o 


2  a 


1 


«=  o 


-  bj/  /0  , 


-b2/  A  . 


Thus,  the  equation  (7.9)  reduces  to 


'  »  1  2 

(7.10)  -j(^-  -  -j-)  sin  2  +  —  sin  *  —  cos  $  . 

A  A 

If  we  further  specialize  and  assume  that  fl  *  C  ,  we  see  from 

(7.10)  that  the  equation  we  now  must  solve  is  merely 

(7.11)  tan  <i  =  b 2 / b ^ 
which  has  solutions 

'l  62 

cos  41  ■  t  —  *  sin  ~ 

where  V  is  the  estimate  of  velocity  defined  in  (7.2).  From  (7.7.4) 
we  see  that  fl  = C  implies  that  the  function  we  must  extremize  is  merely 

and  we  obtain  formally  from  (7.9.1)  the  answer  (which  was  obvious  to 


begin  with) 
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V 


2 


V 


1 

/T 


VIII.  DETERMINATION  OF  A  POSTERIORI  DISTRIBUTION  OF  HEADING  AND 
VELOCITY 

If  the  0.  S.  is  making  observations  in  an  area  which  is  dense 

with  targets,  such  as  a  shipping  lane  or  a  point  of  concentration 

like  the  Straits  of  Gibralter,  and  the  frequencies  of  all  types  of 

ships  and  their  velocities  over  a  recent  period  of  time  are  known 

then  we  can  make  tabulations  of  the  percentages  of  ships  of  various 

types  proceeding  at  different  velocities  in  certain  headings.  This 

characterization  of  the  area  we  assume  can  be  accomplished  with  an 

a  priori  distribution  of  speed  and  heading.  The  a  priori  density 

2 

of  speed  f2(s)  we  take  to  be  N(sQ,d  )  with  the  mean  speed  sQ 
known  as  is  the  standard  deviation  d.  The  a  priori  density  of 
heading  f  is 


fw(h) 


fQ(h  -  w  +  2  ) 
fQ(h  -  w) 


-IT  <  h  <  -IT  +  W 
-H  +  W  <  h  <  TT, 


where 


(8.1) 


fo(h) 


Pjexp 


(^4)  h)-><h<-. 


with  p^,  for  i  *  1,2,  and  the  angle  w  known.  For 

example,  w  is  the  angle  the  shipping  lane  makes  with  the  direction 
of  travel  with  the  0.  S. 
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In  order  to  visualize  the  behavior  of  the  density,  we  present 
a  typical  graph  of  (8.1): 


FIGURE  3 


In  order  for  (8.1)  to  be  a  bona  fide  density,  we  must  have 


TT 

fQ(h)dh 


4. 


By  direct  integration,  this  is  equivalent  to  p^  >  0  and 


%  -  ^  -  >sj  +  p2o2/iv I*  - . 

Following  the  usual  notation,  <t>  denotes  the  cumulative  distribution 
~  function  of  the  standard  N(0,1)  variable  which  again  violates  our 
own  convention  that  Greek  letters  are  vectors. 

Note  that  the  angle  w  can  be  chosen  so  as  to  vary  the  heading 
of  the  greater  percentage  of  ships  away  from  or  toward  the  direction 
of  travel  of  the  0.  S.  Further,  the  constants  p^  can  be  chosen  so 
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as  to  make  the  greater  percentage  of  ships  heading  easterly,  say, 
rather  than  westerly  and  of  course  either  may  be  zero. 


Thus,  having  obtained  an  observation  of  the  random  variable 
(<$,V),  call  it  (4>,v),  we  than  may  obtain  via  Bayes'  theorem 
the  a  posteriori  distribution  of  the  heading  and  speed  (h,s), 
call  it  q,  namely 


(8.2)  q (h,s | $ ,v) 


g($,v|h,s)fw(h)f2(s) 


/•V  2v 

J  J  f (<|>,v|h,s)fw(h)f2(s)dhds 
0  0 


This  integration  can,  in  general,  only  be  carried  out  by  machine 
computation  and,  moreover,  is  suDject  to  all  the  objections  that 
one  can  raise  concerning  the  use  of  prior  distributions.  However, 
in  certain  special  cases  the  density  (8.2)  could  then  be  used  to 
study  the  distribution  of  the  true  heading  and  the  true  speed,  given 
the  observations,  with  meaningful  interpretation  of  the  results. 
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